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Abstract 

Neutrinoless double beta decay can be induced by various mechanisms, involving both light 
and heavy exchange particles such as light Majorana neutrinos, nearly-sterile neutrinos, R-parity 
violating supersymmetric particles, to name a few. It has been pointed out that, typically, the 
contribution due to the exchange of heavy particles between the two nuclei, with masses much 
larger than the typical exchange momentum ~ 100 MeV, is subdominant with respect to the light 
neutrino one. In fact, the former is severely constrained if the same interaction induces light 
neutrino masses. A relevant exception to this generic conclusion is when the contribution to the 
light neutrino masses cancels out. Here, we focus on this case, specifically in the context of seesaw 
models. We perform a general analysis without restricting the study to any particular region of the 
parameter space, although interesting limits associated to inverse and extended seesaw- like models 
are discussed in more detail. It turns out that, once the relevant experimental constraints and 
one- loop corrections to neutrino masses are taken into account, the heavy neutrinos can dominate 
the process only in one of those two limits. For the inverse seesaw we find a very constrained 
allowed region of the parameter space, with heavy neutrino masses around 5 GeV. In general, the 
extended seesaw allows for a larger region, with a very hierarchical heavy spectrum which has 
neutrinos above and below ~ 100 MeV. 
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I. INTRODUCTION 



The existence of neutrino masses, strongly supported by neutrino oscillation experiments, 
is the first experimental evidence of physics beyond the Standard Model (SM). Furthermore, 
the fact that neutrino masses are smaller than the mass of the other SM fermions by several 
orders of magnitude calls for a "natural" New Physics (NP) explanation. Most of the 
models, including the very popular seesaw [IH1] ones, assume that the lepton number is not a 
conserved symmetry and light neutrinos are Majorana particles. An interesting experimental 
window to search for NP gets opened: lepton number violation processes, highly suppressed 
in the SM, among which neutrinoless double beta decay (Oz//3/3 decay) experiments are the 
most promising. In combination with the information coming from neutrino oscillation 
experiments, absolute neutrino mass experiments, precision measurements, and cosmology, 
0u(3(3 decay experiments can give us precious clues in order to identify the mechanism 
responsible for the neutrino mass generation and provide a complementary way to look for 
NP, possibly not otherwise accessible at the LHC. 

Although NP is necessary in order to have Ou/3/3 decay, its effects are indirect as the light 
neutrinos generically dominate the process in most of the models. The key point is the fact 
that the light neutrino contribution and the NP one are usually related through the neutrino 
masses: therefore the NP contribution, suppressed by being short-range, is constrained and 
is generally subdominant [S]. The question if a measurable direct contribution to the Qi//3/3 
decay rate coming from NP is theoretically and phenomenologically viable is thus very 
interesting. This question has been addressed recently in Refs. P, [7] in the context of 
different type I seesaw models. In these publications, a relevant exception to the argument 
above has been pointed out: the case in which the tree level light neutrino contribution, 
induced by NP, partially cancel. In this case, it is found that the direct NP contribution to 
the process is indeed relevant and can be as large as current bounds. However, no detailed 
discussion about the correlation among the light and heavy contributions once the one-loop 
corrections to neutrino masses are included in the analysis is given. The main goal of this 
work is to analyze to what extent this is possible in the general framework of type I seesaw 
models, once the relevant one-loop corrections and experimental constraints are considered. 

We will first very briefly review the aspects of the Ovflfi decay phenomenology relevant 
for our analysis. Considering a general parameterization of the neutrino mass matrix and 
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without restricting the analysis to any region of the parameter space, we will then study 
under what conditions the light neutrino contribution can be cancelled at tree level. We will 
include the one-loop corrections and study if the heavy neutrinos can give a dominant and 
measurable, i.e. at reach in the next-to-next 0^/3/3 decay experiments, contribution to the 
process. We will show that, even when the tree level cancellation takes place, the light and 
heavy contributions are not completely decoupled once the one-loop corrections are included 
in the study and a dominant heavy contribution may occur only in specific regions of the 
parameter space. 

This paper is organized as follows. In section II, we briefly review the Ov/3/3 decay 
phenomenology in the general context of seesaw models, introducing the notation and the 
Nuclear Matrix Elements (NME) we will use. In section III, the parameterization of the 
neutrino mass matrix is presented, distinguishing some relevant limits and their relation 
with well known models, such as the inverse and extended seesaw ones. Section IV is 
devoted to the study of the cancellation condition of the light neutrino contribution and 
its tree level consequences on the heavy neutrino sector. Section V is dedicated to the 
analysis of the relevant corrections: higher order corrections to the seesaw expansion and 
one-loop corrections. The combined analysis of the 0^/3/3 decay phenomenology, when these 
corrections and the relevant experimental constraints are taken into account, is presented in 
section V. Finally, in section VI we draw our conclusions. 

II. DOMINANT HEAVY NEUTRINO CONTRIBUTION TO 0^/3/3 DECAY? 

Generically, in the context of seesaw models, the contribution to neutrinoless double beta 
decay from NP at scales heavier than the exchanged momentum (~ 100 MeV) can not be 
dominant and the light neutrinos typically dominate the process, as shown in Ref. 

Let us very briefly review how the above result is obtained and the possible exceptions. 
Following the notation in Ref. [5j, and restricting the study to type I seesaw models, the 
Ou/3/3 decay rate can be written as 



where Gqi is a well known kinematic factor, U is the unitary matrix which diagonalizes 
the complete neutrino mass matrix both for active and sterile neutrinos, rrij are the corre- 
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(1) 
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sponding eigenvalues, i.e., the neutrino masses, and Ai 0u/313 are the Nuclear Matrix Elements 
(NME) associated to the process. The sum should be made over all the neutrino masses, 
including the heavy ones. 

The NME can be computed using different methods, the main two being the quasi- 
particle random phase approximation (QRPA) [El [9] and the interacting shell model (ISM) 
[Hit ITT] . In this work we will make use of the NME data presented in Ref . [5] and available in 
Ref . [12] . They were computed for different nuclei in the context of the ISM as a function of 
the neutrino mass, something very convenient for our analysis. We use a notation in which 
the dependence on the neutrino propagator is included on A / l 0iy/3/3 (m J '), in contrast with the 
notation usually adopted in the literature where the propagator is expanded to factorize the 
mass dependence. In Fig. 1 of Ref. [5] the NME dependence on the mass of the neutrino 
mediating the process is depicted, showing two different regions separated by the scale of 
the process ~ 100 MeV: 

• Below the 0u/3/3 scale, the NME reach their maximum value and are mainly indepen- 
dent of the neutrino mass. For rm < 100 MeV: M^^rm) = M 0u,3f} (0). 

• The NME corresponding to neutrinos much heavier than 100 MeV are suppressed with 
the heavy neutrino masses and scale as M° v ""(mi) oc 1/mf. 

This behavior of the NME, showing two clearly different regimes, can be easily understood 
expanding the propagator of the neutrino mediating the process. The transition region 
around 100 MeV is well described in Fig. 1 of Ref. [5] since no assumptions have been made 
on the neutrino masses in the NME computation. 

We can distinguish the following two contributions to the Oz//3/3 decay amplitude: 

3 extra 

A oc rrnUlM^imi) + ^ mjU^M ^ '(mi) , (2) 

i=l I 

the first term corresponding to the mostly-active neutrino contribution and the following to 
the extra states of the model. Here and throughout the text, we use capital letters to denote 
the mass index of the mostly-sterile states and lowercase letters for that of the mostly-active 
states. 

On the other hand, since a Majorana mass coupling for the active neutrinos is forbidden by 
the gauge symmetry, the diagonalization of the complete mass matrix leads to the following 
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relation: 

3 extra 

X>E£ + ]£ m/E& = 0. (3) 



i=i i 

This equation, which relates the light and extra degrees of freedom of the model, should be 
always fulfilled at tree level and plays a fundamental role in the phenomenology of 0^/3/3 
decay. 

For extra states with all the masses well above 100 MeV, using the relation given in 
Eq. (|3]), the contribution to Ov/3/3 decay in Eq. (p| can be recast as 

heavy 

A oc m i u ei (M Ou ^{0) - 

I 

heavy 3 
I i=l 

where we have used the fact that Ai 0l, PP(0) ^> A4 0u/3 ^(mi). The contribution from the light 
active neutrinos thus dominates. A similar argument applies to models which implement the 
Type II and Type III seesaw [3] , and more generically to models in which heavy sterile neu- 
trino mixing with v e is introduced. As sterile neutrinos contribute to light neutrino masses, 
Y2i m iUli is constrained by the value of the light neutrino masses while their contribution 
to QvflP decay is suppressed by A4 0ul3 P(mi) making it subdominant, at least if fine tuning 
is not invoked as we will see in the following. 

These considerations apply generically to models with extra sterile neutrinos but there 
are some notable exceptions. 

• The case of extra states below and above 100 MeV. In this case Eq. ^ can be rewritten 

as 

(3 light \ heavy 

i=X I ) I 

3 light \ 

^^ + £m^U^(0), (5) 

. 4=1 / / 

and the new states below 100 MeV may give the dominant contribution. Notice that if 
all the extra states are below the Ou/3/3 scale the cancellation driven by Eq. ^ forbids 
the process. The same behaviour as in this type-I seesaw realization, with extra states 
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below and above the 0vf3f3 scale, can take place in mixed type-I/type-II (type-I/type- 
III) scenario. In these scenarios NP above the Qvflfl scale, such heavy sterile neutrinos 
or heavy triplets in the just mentioned cases, is needed to avoid the cancellation, while 
the extra "light" sterile neutrinos may give an important contribution to neutrino less 
double beta decay [5]. 

• Additional contributions to neutrino masses. In this case the mass relation becomes 

3 heavy 

mU^i + Y m I U eI = m LL, (6) 
i=l I 

where tull is an effective Majorana mass term generated for the active neutrinos by 
some other mechanism, mn and ^^ eavy miU^j could be very large and cancel nearly 
exactly, keeping light neutrino masses under control. In this way, even with the con- 
tribution to Ov/3/3 of the heavy states being weighted by the corresponding NME, a 
dominant effect could arise. However, it would have to overcome the suppression com- 
ing from the NME (A / t 0iy/3/3 (0)/A^ 0y/3/3 (m / ) > 1) and a very high level of cancellation 
among ^^ eavy mjU^ and ttill in Eq. would be required. This possibly implies an 
uncomfortably high level of fine tuning and will not be studied in this work. 

• A cancellation in the light neutrino contribution: Y2i=i m i^ei = 0- If this cancellation 
takes place, the heavy neutrinos would trivially dominate the process (at least, at tree 
level) . 

In this work we are going to focus on this last possibility. This relevant exception was 
studied in Refs. [HI El H3] and not contemplated in Ref. [5J. Of course, this cancellation 
in the light contribution could be obtained invoking some symmetry, and the simplest one 
in this context is the lepton number. The well known inverse [H] or linear [15] seesaw 
models, that involve small violations of the lepton number, could in principle implement 
this scenario. However, generating a measurable heavy neutrino contribution to the Qi//3/3 
decay is not trivial even in these models. First of all, Ov/3/3 decay is a lepton number 
violating process and consequently is expected to be suppressed in this context. Moreover, 
the suppression of the heavy neutrino contribution with the NME (~ 1/mj) makes having 
very low scale heavy masses unavoidable in order to obtain a relevant effect. This possibility 
has been recently explored claiming that the heavy neutrinos could be very relevant for some 
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particular neutrino mass textures [6j [7]. Indeed, the main goal of this note is to check to 
what extent this is possible, paying special attention to the stability of the light neutrino 
masses under one-loop corrections and their contribution to the Qv/3/3 decay. 



III. THE MODELS 

We will focus on the study of SM extensions which consist in the addition of n+n' fermion 
gauge singlets, N^, to the SM particle content without imposing lepton number conservation, 
whose Lagrangian is 
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JSf = S? m + ^kin - ^N l M ij N'r - {Y) ia N^L a + h.c. 



(7) 



where j£?sm is the SM Lagrangian and =S?km ar e the kinetic terms of the new fields iVj. Here, 
and in the rest of the paper, the sub index a denotes flavour (a = e, /i,r). Without loss of 
generality, the neutrino mass matrix can be expressed as 



( eY 2 T v/V2 \ 



-1 


< 


m T D \ 




^ m D 


M J 



•U.= Y lV /V2 fi' A = | " | • (8) 

\eY 2 v/V2 A T ii j 

Here Y\ and Y 2 are the n x 3 and n'x3 matrices that form the Dirac block mp. The Majorana 
submatrix M is composed by //, and A: n x n, n' x n' and n x n' matrices respectively. 
Another helpful, and widely used, basis is the one in which the Majorana submatrix for the 
sterile neutrinos is diagonal, which we will denote with "tilde" in the following discussion. 
In order to illustrate the relation between these two basis, let us consider the n — n' — 1 
case. The Majorana submatrix M can be diagonalised as 



OM v O T = M v , 



O 




(9) 



being A a 2 x 2 orthogonal matrix with rotation angle 1 

fi'-H+ V4A 2 + (fx' - /i) 2 



tan 9 



2A 



(10) 



1 For simplicity, all the Majorana submatrix parameters in M have been considered real. 
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The Majorana masses M\ and M 2 and the Yukawa couplings Yi are then given by 

M 2 ,i = \ (V + fi ± V4A 2 + (//' - /i) 2 ) , 
fi = Yi cos - eY" 2 sin , 

F 2 = Yisin0 + eY 2 cos0 . (11) 

Of course, the analysis can be performed in any basis, but we will mainly work in the one 
in which the neutrino mass matrix is given by Eq. (|8]). 

Notice that the mass matrix given in Eq. (JsJ) is completely general. A particularly inter- 
esting set of models, included in Eq. (J8J), are those studied and summarized in Ref. [16] and 
which include the so-called inverse or multiple seesaw models [T ^ IT7HT§] . Lepton number 
is assumed to be a good global symmetry only broken in the neutrino sector through the 
small lepton number violation terms e and/or /i and/or //. 

The light masses are proportional to e and fi as long as Mj ^> mr>. In this way, their 
smallness could be due to the one of e and fi, small parameters which break the lepton 
symmetry. Therefore, the scale of the NP, given by A, can be lowered to the TeV scale or 
even below. This allows sizable NP effects coming from the dimension six operator which, 
contrary to the dimension 5 one, does not present an extra suppression with e and fi as 
it does not violate lepton number [16] . Lepton conserving processes very supressed in the 
SM as the rare decays are very promising channels to probe this kind of NP. Interesting 
recent analysis of \i — > ej, fi — >■ eee or fi — > e conversion in the context of low scale small 
lepton number violating models can be found in the literature [201422] . In Ref. [6] these sort 
of models are studied in the context of the 0^/3/3 decay using a different parameterisation 
based on the Casas-Ibarra one [23J, which parameterises the neutrino mass matrix in terms 
of the light and heavy masses, the U a i matrix and an orthogonal matrix R. Notice that the 
parameterisation considered here is totally general and includes the Casas-Ibarra limit in 
which an approximate decoupling of light and heavy sectors is assumed. 

In any case, we will not restrict our study to any particular value of the parameters or, 
in other words, to any of the above mentioned specific limits. Nevertheless, for simplicity, 
we will consider the case in which only two fermion singlets (n = n' = 1) are added. The 
conclusions obtained in this work can be extended to models with larger number of right 
handed neutrinos, barring fine tuned cancellations. On the other hand, from neutrino os- 
cillations we know that it is not easy to accommodate the experimental data in the region 
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of the parameter space between the limits: M, ^> % (seesaw limit) and Mj <C % (pseu- 
doDirac limit). In fact, in Ref. |24J it is shown how the constraints from neutrino oscillation 
experiments leave those limits as the only allowed regions for n — n' — 1 and M\ = M 2 . 
The region of the parameter space in between is ruled out and only the pseudoDirac and 
seesaw limits survive. Reasonably extrapolating these results to the more general case with 
Mi 7^ M 2 studied here, leaves the seesaw limit (Mj 3> rfi£>) as the only relevant part of the 
parameter space in the Oz//3/3 decay context 2 . From now on, we will focus on the seesaw 
limit. Notice, however, that this does not necessarily mean that Mj have to be at the GUT 
or the TeV scale and can be considerably lighter [25T - T27] . 



IV. LIGHT NEUTRINO MASSES AND Ov/3/3 DECAY 

For Mj ^> rfiD, the light neutrino mass matrix is given at tree level by 

m tree ~ -m T D M- x m D ~ ^f—— (jiY^Y x + e 2 // Y 2 T Y 2 - Ke{Y?Y 1 + Y?Y 2 )) , (12) 

where and M are the 2x3 Dirac and 2x2 Majorana sub-matrices respectively in 
Eq. (|sj) for n = n' = 1. Here, we have performed the standard "seesaw" mn/M expansion 
keeping the leading order terms. We will discuss later if the higher order corrections can be 
relevant. The contribution of the light mostly-active neutrinos to the Oz//3/3 decay amplitude 
is proportional to the "ee" element of this effective mass matrix as 

A ligM oc J2 m * U eiM°^(0) « - KM- 1 m D ) ee ^°^(0) = 
i=i 

Therefore, the light neutrino contribution is strictly cancelled as long as the parameters of 
the model satisfy the following relation 

HYI + tY 2e (efi'Y 2e - 2AY le ) = . (14) 

This condition is fulfilled for 

e = ^ = 0. (15) 

2 Of course, the Dirac limit will not be considered in this analysis where the Ov/3/3 decay phenomenology is 
studied. 
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Of course, it may also be satisfied for other choices of parameters, but e = /x = is the 
most stable one under radiative corrections and higher order terms in the expansion, as we 
will show later. From now on we will assume that this cancellation condition is fulfilled. 
Obviously, setting e and \x to zero leads to vanishing tree level active neutrino masses as 
well. However, the light neutrino masses can be generated at one-loop as we will see. 

One could naively think that taking into account Eq. (|3]) would lead us to the same 
cancellation for the heavy neutrinos (see Eq. (|2j), however, the dependence of the NME on 
mj avoids a complete cancellation, if the heavy neutrinos are not very degenerate. 

When the heavy neutrinos are above the 0i>(3(3 scale, 777,4, m$ ^ 100 MeV, the heavy 
contribution to the 0v/3f3 decay amplitude can be approximated as 

extra 

A extra oc rnjU^M^imj) oc - (m T D M^m D ) ee (16) 



1 



,2 



(/x 3 + A 2 (2 M + mQ) Yl - 2eA (A 2 + /Z 2 + p 2 + /V) Y le Y 2e + (// 2 + A 2 (p + 2fJ)) e 2 Y 2 e 

2(A 2 - W ') 3 

which reduces to 

Aextra OC . (17) 

if the light neutrino contribution is cancelled (e = fi = 0). Apparently, the above expression 
indicates that for large values of // and/or small enough A the heavy neutrinos may give a 
relevant contribution to the Ov/3/3 decay at tree level. At this point two interesting limits of 
Eq. pi) arise: 



Extended seesaw limit (ESS limit): pi! ^> A, mp. In view of Eq. (17), this possibility 
appears quite appealing. This limit matches the so-called extended seesaw models [2H] 
and corresponds to a hierarchical spectrum for the heavy neutrinos: 

m 4 w Mi w -A 2 //i', U ei w Y le v/V2A, q ^ 

(loj 

m 5 » M 2 W /i', «7 e5 w Y u v/y/2fi', 

where we also show the corresponding mixing with the active neutrinos. In this regime, 
the lightest of the two heavy neutrinos dominates the heavy contribution. Moreover, 
for large enough values of //, 7774 becomes lighter than 100 MeV, the NME takes its 
maximum value and the heavy contribution to the Ov/3/3 decay becomes independent 
of A: 

A extra oc U 2 4 m 4 M°^(0) « -|^M^(0) . (19) 
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• Inverse seesaw limit (ISS limit): A 3> //, mo. This limit corresponds to one of the 
Minimal Flavour Violation models (MFV) studied in Ref. [16J. It is also related to the 
case analyzed in Ref. [5], where a different parameterization is used. In this case the 
heavy neutrino spectrum is quasi-degenerate, forming a quasi-Dirac pair: 

m 4 ps — m 5 ps Mi ps — M 2 ps A, [7 e4 ps [7 e5 ps Y" le w/2A, ^ 
AM = |M 2 | — |Mx| ps //, 

and we can expect lepton number violating processes such as neutrino less double beta 
decay to be controlled by //. 

If all the heavy neutrinos are located below the 0is(3/3 scale, a cancellation driven by 
Eq. ^ is expected at tree level, as we have already mentioned. This cancellation applies in 
general as long as all the heavy neutrinos are in the light regime, including the two limits 
distinguished above. 



The approximation made in Eq. (16), M 0ul3 ^(mj) oc l/mf, does not apply if one of the 
heavy neutrinos (or both) is lighter than (or close to) ~ 100 MeV. However, as we have 
already commented, we will not restrict the analysis to any particular value of the sterile 
neutrino masses. This is the reason why we have made use of a numerical computation for 
the NME in which no approximation for the neutrino mass dependence has been considered. 
Notice, for instance, that the phenomenology for heavy masses around 100 MeV can be very 
interesting and the approximation A4 0lJ ^^(mi) oc 1/mj is not very accurate in that region. 

In summary, at tree level the light neutrino masses are independent of // (and A) for 
e = /j, = 0, being actually zero. However, lepton number violation processes such as Qi/fiP 
decay are sensitive to these parameters and \J in particular. The idea behind Ref. [BJ [7] is to 
exploit this apparent decoupling between the heavy and light contributions in order to have 
a measurable effect in the 0u/3/3 decay coming from the heavy side. In the following, we will 
check if a heavy dominant contribution is really possible once the relevant corrections and 
experimental constraints are taken into account. 

V. HIGHER ORDER CORRECTIONS IN THE SEESAW EXPANSION 

Only the leading order in mo/M has been considered in the expansion performed in 



Eq. (12). We now check if the higher order corrections may induce any relevant effects 
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once the tree level cancellation for the light masses takes place. The next to leading order 
contributions to the light neutrino masses can be written as [29J: 



§m = \n tree m} D M 2 m D + ^ (m tTee m ] D M 2 m D ^j , (21) 

where m tree is the leading order contribution given by m tree = —nv[ ) M~ 1 m£). As they are 
proportional to the leading order active neutrino mass mtree, they are completely irrelevant 
for fi — e — 0. In fact, the light neutrino masses vanish for fi = e = at all orders in the 
expansion [2U| I3U] . Contrary to the \i — e — case, other choices of the parameters which 



satisfy the cancellation condition given in Eq. (14) are flavour dependent, giving as a result 
non-vanishing higher order corrections. 

On the other hand, the factor m} D M~ 2 mD/2 is nothing but the coefficient of the effective 
d = 6 operator obtained when the heavy neutrinos are integrated out of the theory |31j. 
This coefficient, which induces deviations from unitarity of the 3x3 lepton mixing matrix, 
is independent of // when the light neutrino cancellation (/i = e = 0) takes place. Therefore, 
for fM — e — 0, the d = 6 operator does not introduce any relevant //-dependent deviation 
from unitarity and \j! can escape from the corresponding constraints [32, 33J, even if fi' ^> A. 



VI. ONE-LOOP CORRECTIONS 



The one-loop corrections can be of two different types: renormalizable (i.e. the running of 
the parameters) or finite. In this section we will study both, starting with the renormalizable 
corrections. The analysis will be done after electroweak symmetry breaking (EWSB). 



A. Renormalizable one-loop corrections 

We are mainly interested in the running behavior of the parameters /i and e, since they 
drive the light neutrino mass cancellation, and /i' and A which are the key parameters 
associated with the heavy contribution. We have performed the computation in the basis 
in which the neutrino mass matrix takes the form given by Eq. (|8]), in such a way that the 
one-loop running equations [331436] for these parameters can be directly obtained: 
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Q d(eY 2a 



T~\g 2 - \g' 2 \ Y 2a - I Y w (jyfYfifr - Y* p Y la ) + \e 2 Y 20 Y* p Y 2a 



dQ (4tt) 2 

Q % = J^yV Y i^ + ^Y*,Y w ], 

Q % = jhf [A Y ^ + e ^' Y ^ Yw + * Y; ^ + A e Y; ^ ' (22) 

with T = Tr (3Y^Y U + W ] d Y d + Y?Y t + Y^Y^j and being g and g' the 5Z7(2) L and U(1) Y 
gauge coupling constants of the SM. We do not need to solve the equations to realize that 
the effect of the one-loop renormalizable corrections to \x and e is suppressed by the tree 
level values of e or //. This means that the cancellation of the light active neutrino masses is 
stable under one-loop renormalizable corrections, as expected as a Majorana mass coupling 
for the active neutrinos is not allowed at tree level. For vanishing e and // at tree level, the 
light neutrino masses keep being zero independently of the running of the parameters (even 
for huge tree level inputs of //). This is no longer true once the finite corrections are taken 
into account as we show in the next subsection. 



B. Finite one-loop corrections 

Indeed, after EWSB, a Majorana mass for the active neutrinos is generated through 
finite one-loop corrections. Of course, the other Yukawa and Majorana couplings among 
the active and sterile neutrinos also get finite corrections, but their contribution to the light 
neutrino masses vanish for \i — e = 0. Therefore, the dominant contribution to the light 
neutrino masses comes from the Majorana mass generated for the active neutrinos and is 
given by |37H3H] 

1 T f 3 In (M 2 /M|) In (M 2 /M 2 H )\ 

where and M are the Dirac and Majorana sub-matrices respectively, Mz is the mass of 
the Z boson and Mh the Higgs boson mass. Notice that no expansion has been performed 
in order to obtain this result. The structure of the correction is similar to the tree level 
masses but in this case no cancellation takes place for \x = e = 0. 
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FIG. 1: The coloured region in the //-A plane corresponds to \5rriLL (fi = e = 0) | > 0.1 eV. 
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respectively. 



Notice that the above expression is valid in any basis. In particular, Eq. (23) can be 
conveniently written in the [i = e = limit as 



5m 



LL 



1 Y{Yi 



~ r 2 



3Mi In Mi /M 



Mi In [Mi/M 2 H 



Mi/M 2 z - 1 



+ 



+ 



(4tt) 2 2 

3M 2 In f M|/Af| ) M 2 In ( M|/M| 



M|/M| 



+ 



M|/M| - 1 



Mi/M 2 H - 1 



sin 2 fl 



cos 2 ^ 



(24) 



where the M2,i are the eigenvalues of the Majorana mass term given by Eq. (11), and 9 



is the rotation angle given by Eq. (10), both evaluated for e = fi = 0. In Fig. [T]we show 
the region of the parameter space //-A given by \5rriLL (/x = e = 0) | > 0.1 eV for different 



values of the Yukawa couplings. In order to understand better the implications of Eq. (24), 
we have obtained approximate expressions for two relevant limits: 
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A > n', M H , M z . We have 



5m 



LL 



1 Y{ Y 1 Mfj + 3M| 



(4tt) 



A 2 



(25) 



As we have already discussed in Sec. [TV] this case is included in the ISS limit and 
corresponds to a MFV model in which /x, e and // are lepton number violation pa- 
rameters. What we observe here is that, although the tree level light neutrino masses 
cancel for e = /x = 0, they are generated at one loop and are proportional to the 
only lepton number violation parameter different from zero, /x', as expected since the 
neutrino masses also violate this symmetry. 



/x' ^> A 3> Mh, Mz- In this case one finds 



5m 



LL 



1 Y?Yi 



(47T) 



3Mj 



In 



A*_ 
Ml 



+ 



M 



xx 



In 



A 4 
Ml 



(26) 



This case is included in the ESS limit discussed in Sec. |IV| Here, the one-loop light 
neutrino masses depend mildly on A and are suppressed by /x'. Again this can be 
understood in terms of a lepton symmetry: /x' is suppressing the violation of Lepton 
number at low energies in such a way that in the limit /x' — > oo the symmetry is 
completely restored in the effective theory. 



In the next section we will study the phenomenological consequences of Eq. (24) in the 



context of the Ov/3/3 decay without considering any expansion on the parameters. It is 
important to remark here that, once the tree level cancellation takes place, only one mass is 
generated at one-loop and at least two light masses are necessary to explain the light neutrino 
spectrum obtained in neutrino oscillation experiments. This is easy to solve: simply adding 
another fermion singlet to the model would allow to generate the necessary extra light mass. 
Although, for simplicity, we will keep studying the simpler case with only two extra sterile 
neutrinos, the conclusions extracted from the analysis can be extended to models with more 
than two singlets, barring fine tuned cancellations. 
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VII. NP DOMINANT CONTRIBUTION TO Oi/fiP DECAY AND ONE-LOOP 
NEUTRINO MASSES 



Once the relevant one- loop corrections are taken into account, the Lagrangian is modified 

to 

if = ^sm + ^kin - \MM l3 Nc - {5m LL ) aP P^v^ L - {Y) ia N$L a + h.c. . (27) 

Consequently, Eq. pj), which comes from the diagonalization of the neutrino mass matrix, 
is also modified to the following one-loop version: 

light extra 

^iUei + miU *I = ( 6m LL)ee- (28) 

i I 

In the case of interest, when e = fi — 0, the light neutrino masses associated to the mostly- 
active neutrinos are determined by 5m ll- The tree level condition YIT^^ m i^ei = remains 
true at one-loop level but, as discussed in Sec. [TV] heavy neutrinos could have a sizable effect 
to 0^/3/3 decay thanks to the NME dependence on the heavy masses. However, at one-loop, 
the NP contribution to the OvftP decay and the light neutrino masses are related as they 
depend on the same parameters, in particular //, A and Y\ a . Their decoupling achieved 
at tree level does not remain true once radiative corrections are included. Consequently, 
the heavy parameters can not be chosen arbitrarily as to dominate Oz//3/3 decay but are 
constrained by light neutrino masses, which also contribute to the process. 

In principle the radiative corrections dependent on the transfered momentum p have also 
to be considered. This corrections are of two types: i) proportional to p or p 3 ; ii) dependent 
on p 2 . The first come from the W and charged Golstone boson corrections to the neutrino 
propagator and vanish after performing the integration over p in the Ou/3/3 amplitude. The 
second are associated to the Z and Higgs boson corrections to the propagator and are 
negligible in the region of heavy masses under consideration. 

In the rest of the section, we will study under which conditions it may (or may not) be 
possible to have a dominant heavy neutrino contribution. We will pay special attention to 
the impact of the one-loop corrections and the experimental constraints on the parameters 
of the model. To illuminate the interplay among all these factors, we will first analyze the 
particular case of Y\ a = 10~ 3 showing our results in Fig. [2j 

First of all, we are assuming that the model under consideration provides the dominant 
source of light neutrino masses. In principle, they should be in agreement with neutrino 

16 



> 

CD 



&0 

o 




5 6 7 

Log [ A/Me V] 



FIG. 2: Impact of one-loop corrections to decay for Y\ a = 10 -3 . The red band is the 

95%C.L. allowed region for the one-loop generated light neutrino masses bounded by cosmology 
and neutrino oscillations. The orange band is the 95%C.L. region of the parameter space in which 
the heavy neutrino contribution is between the present bound from EXO and the sensitivity of the 
next to next generation of experiments. Blue (green) stands for to region in which the ratio r 
between the heavy and light contributions is r > 5 (1 < r < 5). The grey region inside the dashed 
black line is the parameter space ruled out at the 95% C.L. by the constraints on the mixing. 

oscillations data but, since we are generating at one loop only one light neutrino mass we 



only impose a conservative lower bound on the non trivial eigenvalue of Eq. (24) given by 



the solar splitting, y'Am^. Moreover, the absolute neutrino mass scale experiments impose 
an upper bound on the same combination of parameters and as reference value we take the 
95%C.L. upper bound on the light neutrino masses from cosmology [41j . m u = 0.58 eV. Since 
we are analyzing the case in which the tree level active neutrino masses cancel (e = /x = 0), 
these bounds can be directly translated into bounds on // and A as a function of Y\ a . They 
are shown in Fig. [2] as the red band. The Higgs mass, m#, has been fixed to 125 GeV 
in all the calculations, as suggested by the recent LHC results (US IHj. Notice that if no 
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lower bound is imposed, as it would be the case if the light neutrino masses come from some 
other mechanism, the outer region of the red band would not be excluded. Our conclusions 
remain valid also in this case, as we will discuss later. The constraint on the light neutrino 
masses shown in Fig. [2] can be understood analytically taking into account the approximate 
expressions derived in the previous section. In the ISS limit a/ Am^ < Stull < 0.58 eV 



scales as /i'/A 2 in agreement with Eq. (25). For // 3> A, in the ESS limit, it becomes mainly 



independent of A according to Eq. (26). 

The heavy neutrino contribution to Ou/3/3 decay, given by -Aheavy CX- 
J2i=4 5^ei m i-^° u ^ l3 ( m i)^ can be computed diagonalizing the mass matrix in Eq. (JsJ) 
and using the NME data calculated as a function of the neutrino masses [12]. The 
diagonalization can be easily performed in the e = [i = limit. This contribution has to 
respect the present experimental bound and, in order to be phenomenologically interesting, 
should be at reach of future 0^(3 (3 decay experiments such as CUORE [45J, EXO |46j, 
GERDA [Sj, KamLAND-Zen gg], MAJORANA 03], NEXT |5(l or Super-NEMO [EE]. 
This constraint is shown as the orange band in Fig. § the 95%C.L. region of the parameter 
space in which the heavy neutrino contribution is between the present bound from EXO [46] 
(Oz//3/3 decay in 136 Xe), which using the corresponding shell model NME is \mpp\ < 0.53 
eV, and the future (optimistic) sensitivity of the next-to-next generation of 0v(3(3 decay 
experiments, taken to be mpp = 10~ 2 eV. The shape of the heavy contribution contour can 



also be easily understood from the discussion in Sec. IV The heavy contribution scales as 



yti'/A 4 following closely Eq. (17) until, for // ^> A, becomes independent of A in agreement 



with Eq. (19). In the ISS region, both heavy neutrinos have masses lager than the Qvf3f3 



scale and Eq. (17) holds. As expected comparing Eq. (17) and Eq. (25), in this region 



the slope of the heavy contribution contour is twice the a/Attt^ < Stull < 0.58 eV one. 
For > A we enter the ESS limit and eventually the lightest of the two heavy masses 
becomes lighter than the Ovfifi scale, while the heaviest one is too heavy to give a relevant 
contribution to the process. The "heavy" contribution is dominated thus by the sterile 
neutrino lighter than 100 MeV, for which the corresponding NME takes the maximum 



value, and is independent of A, see Eq. (19). 

Fig. [2] also highlights the region of the parameter space for which the ratio r between the 
heavy and mostly-active contribution to Ov/3/3 decay, defined as r = \A heavy /An ght \, is be- 
tween 1 and 5 (green region) or larger than 5 (blue region). The active contribution is deter- 
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mined by the one-loop correction to the light neutrino masses: A actiive oc (Sttill) ee A4 0u/313 (0) . 
From Eqs. (17), (19) and (24), it is clear that r = \A heavy /An ght \ should be basically inde- 



pendent of the Yukawa couplings. 






5 6 7 

Log[A/MeV] 



5 6 7 
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FIG. 3: Same as in Fig. |2| for Y la = 10" 4 (left), Y la = 10" 5 (center) and Y la = 3 • 1CT 6 (right). 

Finally, the information coming from the experiments that constrain the mixing between 
the active and heavy neutrinos is also included in Fig. [2] The grey region inside the dashed 
line is excluded at the 95% CL. by the constraints on the mixing extracted from weak decays 
and summarized in [52J and non-unitarity bounds [32J E3] • 

As shown in Fig. |2j it is possible to have a dominant and measurable contribution from 
the heavy neutrinos to 0u/3j3 decay keeping the light neutrino masses under control. Ignoring 
for the sake of the discussion the constraints on the heavy mixing, this takes place in the two 
intersections among the red, the orange and the blue regions, which lie in two interesting 



limits already discussed in Sees. Ill and IV 



• i) ISS limit: A ^> //, Y\ a v. The heavy neutrinos are quasi-degenerate, and their contri- 
bution to the process proportional to the splitting, given by //. Once the constraints 
on the mixing, f/ e4 ~ c7 e5 ~ Yi e t>/2A, are properly taken into account the ISS limit is 
ruled out. 

• ii) ESS limit: // ^> A, Y\ a v. In this case, the lightest of the extra neutrinos has mass 
lower than the neutrino less double beta decay exchange momentum and dominates 
the process. 

We have chosen 10~ 3 as the input value of Y\ a in Fig. [2] as an example that allowed us to 
illuminate the discussion. The results for Y\ a = 10~ 2 — 1CF 3 are pretty similar, but we have 
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FIG. 4: Region of the parameter space, M 2 ~ 7715 vs Mi ~ 7714, in which a dominant and mea- 
surable contribution of the heavy neutrinos, respecting bounds from neutrino oscillations, absolute 
neutrino mass scale experiments and weak decays, is feasible. Blue, cyan, green, yellow and red 
stands for Y\ a = 10 -2 , 1CP 3 , 10~ 4 , 1CP 5 and 3 • 1CP 6 respectively. The black lines correspond to 
Mi = 100 MeV and M 2 = 100 MeV. 

checked that for values of the Yukawa couplings larger than 10 -2 , a dominant contribution 
from the heavy neutrinos is not possible and can be at most of the same order as the one 
from light neutrinos. In Fig. [3] we show the plots analogous to Fig. [2] but for smaller values 
of the Yukawa couplings: 10~ 4 (left), 1CT 5 (center) and 3 • 10~ 6 (right). For these values, 
the heavy neutrino mixing is small enough to satisfy the bounds coming from weak decays. 
We observe that the ratio between the light and heavy contributions is independent of the 
Yukawa couplings as expected. However, each of them separately depends strongly on that 
input. The region of the parameter space in which we have a measurable heavy contribution 
decreases with the Yukawa coupling, as it is also the case for the red region, in which the 
light neutrino masses keep being under control. We have checked that, between 10~ 6 and 
10 -8 a dominant and measurable contribution of the heavy neutrinos may still be possible, 
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but the light neutrino masses generated at one-loop are smaller than y'Am^. For values 
of the Yukawa couplings smaller than 10~ 8 the heavy contribution is too suppressed to be 
experimentally accessible. 

The information given in Figs. [2] and [3] is summarized in Fig. [4j where we show the 
region of the parameter space in which a dominant and measurable contribution of the 
heavy neutrinos is possible respecting at the same time the bounds on heavy mixing from 
weak decays and non-unitarity [32j |33j [52] and keeping light neutrino masses in the region 
between ^/Am^ and their upper bound extracted from Ref. [H] . Although the tree level 
cancellation for the light neutrino masses is taking place, once the one-loop corrections are 
taken into account, a dominant contribution from the heavy neutrinos can not occur for 
larger or smaller values of the Yukawa couplings than the ones shown in Fig. |4| Notice 
that this dominant contribution is mainly possible only in the hierarchical seesaw scenario 
mentioned above (|Mi| < 100 MeV <C \M 2 1), where the lightest sterile neutrino gets a mass 
smaller than (or around) 100 MeV and dominates the process. Indeed, this result is not 
surprising: in Ref. [5] it was shown that, in the case in which the cancellation of the light 
neutrinos contribution does not occur, a hierarchical heavy spectrum like this is necessary 
in order to have a relevant contribution from the heavy neutrinos at tree level. We have 
checked in this work that this conclusion, obtained at tree level, can be extended to the 
case in which a cancellation of the light contribution takes place at tree level once the one- 
loop level corrections are included in the analysis. Nevertheless, there is an exception to 
these conclusions. For Y" lQ « 10~ 4 — 10 -5 , there is still a tiny region in which the heavy 
contribution could dominate when the heavy neutrinos are quasi-degenerate and around 5 
GeV (ISS region). 

Finally, a comment is in order: the qualitative conclusions just depicted above are not 
affected if the lower bound on the one-loop light neutrino masses imposed here (Sitlll > 
^Am^) is not assumed in the analysis. In such a case, the allowed regions in Fig. ^become 
a bit larger (vertically) and a dominant heavy neutrino contribution would be still possible 
for Y\ a = 10~ 6 — 10~ 8 . Also in this case a hierarchical spectrum with \M\\ < 100 MeV 
<C IM2I is required with the possible exception of having a quasi-degenerate spectrum with 
I Mi I ~ IM2I ~ 5 GeV (in the same tiny region of the parameter space). This can be easily 
understood from Figs. |2](3j eliminating the lower bound on Stull would mean that the outer 
region of the red bands would not be forbidden any more. 
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VIII. SUMMARY AND CONCLUSIONS 



The possibility of having a dominant contribution from "heavy" neutrinos to Ois/3/3 decay, 
when a cancellation of the tree level light neutrino contribution takes place, has recently re- 
ceived much attention P, [7] . In this work we have carefully analyzed this possibility in the 
general framework of type-I seesaw models. We have considered a general parameterization 
of the neutrino mass matrix which allowed us to explore the whole parameter space, identify- 
ing particularly interesting limits such as the inverse and extended seesaw models. We have 
shown what conditions have to be satisfied for a stable cancellation of the tree level light 
neutrino contribution, allowing the heavy neutrinos to dominate the process at tree level. 
We have studied the relevant corrections that may arise in this context. The finite one-loop 
corrections to the light neutrino masses turn out to be very relevant. Although logarithmic, 
their contribution to the 0^/3/3 decay rate tends to dominate very easily. We have found that 
the heavy neutrinos can give the main contribution to the process only for a very hierarchical 
heavy neutrino spectrum with masses below and above the Ou/3/3 scale ~ 100 MeV, which 
would match an extended seesaw like model. The heavy neutrino contribution is in fact 
completely dominated by the lightest sterile neutrinos with mass < 100 MeV, which is not 
suppressed by the NME. This result coincides with the general conclusions of the tree level 
analysis performed when no cancellation takes place [5]. Quantitatively, we have obtained 
that values of the Yukawa couplings between 10 -2 and 10~ 6 (10~ 8 ) are necessary, if a lower 
bound on the one-loop neutrino masses of a/ 3 m 1 i (no lower bound) is imposed in the anal- 
ysis. We qualitatively agree with part of the conclusions drawn in Ref. [Zj: the extended 
seesaw scenario might accommodate a relevant "heavy" neutrino contribution. Nevertheless, 
our general conclusions clarify an important detail: in Ref. [7j it was hypothesised that this 
may generically happen with all the heavy neutrinos above the Ou/3/3 decay scale while we 
conclude that the heavy spectrum needs to contain states in both regimes, below (or close 
to) and above 100 MeV. 

An interesting exception arises for quasi-degenerate heavy neutrinos with masses around 
5 GeV which may give the dominant contribution in a tiny region of the parameter space for 
Yukawa couplings in the range 10 -4 — 10~ 5 . In agreement with Ref. [6j [7j, we confirm that 
a relevant contribution to the Ou/3/3 decay may come from a seesaw scenario with a quasi- 
degenerate heavy neutrino spectrum, which corresponds to an inverse seesaw like model. 
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However, we also show that this possibility is rather unlikely since it can only take place in 
a very particular and small region of the parameter space. 

We should remark that our analysis was performed considering just two extra fermion 
singlets in the model. But the conclusions obtained can be safely extended to scenarios with 
more fermion singlets barring (more) fine tuned cancellations. 

Note 

During the completion of this work an analysis which considered the generation of 
neutrino masses at the loop-level in inverse seesaw models was presented in Ref . [53] . 
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